Numerical Test of the Energy Scale of Columnar Dimerization in High-Tc Cuperate 

Superconductors 
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Spontaneous bond dimerization (BD) was cojectured to be central to the physics of the cuperate 
superconductors. Based on ladder models and numerical simulations, we examine the quantitative 
issue of whether the energy scales involved are sufficient to influence the mechanisms of pairing, stripe 
formation, pseudo-gap, and quantum critical 'normal' state phenomena. We present preliminary 
evidence that maximum gap scales derived from BD are always too small. 



A mechanism of spontaneous bond dimerization (BD) 
was proposed P, U S 0] as an underlying core phe- 
nomenon in the cuperate high temperature supercon- 
ducting materials. Columnar BD was proposed as the 
mechanism for the formation of stripes in underdoped cu- 
perates (see Fig.l). A quantum critical point associated 
with BD is proposed to be responsible for certain anoma- 
lous thermodynamic and magnetic properties observed in 
the temperature-doping phase diagram of the cuperates. 
Elaborate consequences for impurity and vortex states 
were derived 0, Moreover, the BD is characterized 
by a spin gap and associated confinement of spin-i exci- 
tations, and hence we note that by itself it also provides 
a potential mechanism of pairing. 

The analysis in Refs. Q, H, |^ focussed primarily on 
symmetry aspects of the BCS and BD order parameters. 
Here we focus on quantitative issues. In the cuperates, 
the maximum spin-gap in the superconducting state is 
As = 2 A (0) « 800K ~ J/2 [3, where J is the antifer- 
romagnetic exchange interaction. Pseudo-gap phenom- 
ena persist to temperatures of order 300K. Some evi- 
dence of local stripe order (incommensurate magnetism) 
persists up to temperatures of order 70K in some uir- 
derdoped materials, where spectral signatures of incom- 
mensurate magnetism are seen up to energies of order 
AOmeV (~ @|. Thus, for these phenomena to be 

associated with BD order or fluctuations, this order must 
generate an energy scale, say a spin gap (or spin pseudo 
gap), which is sufficiently large. We will adopt the le- 
nient criterion that this energy scale must be at least of 
order A > 250K « J/6. 

Our main observation is that this aspect of the BD con- 
jecture can be tested by numerical simulations on wide 
ladders (e.g., 6- leg or 8-leg). Indeed, existing numer- 
ical experiments on 4-leg ladders and on small square 
Heisenberg arrays will be shown to already suggest rather 
strongly that the energy scales derived directly from BD 
are coirsiderably smaller than 250K. Pending sharper 
similar results in 6-leg ladders, we doubt the viability 
of BD as a candidate theory of pairing or stripes in the 
cuperate high- Tc superconductors. 

We stress that our analysis bares strictly on the mech- 
anism of spontaneous BD. Associated results of the work 
by Sachdev and collaborators may be relevant to more 
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FIG. 1: Conjectured phase diagram by Sachdev et al. 



general models of stripes or ladders, and hence their rel- 
evance to the cuperate superconductors may hitherto not 
be tied to the mechanism of spontaneous BD. 



REVIEW OF THE COLUMNAR DIMERIZATION 
APPROACH TO HIGH-TC 

Columnar BD was introduced as a bond ordered state 
in 2D Heisenberg models with frustrated interactions 
4]. The caironical example is the model with antifer- 
romagnetic irearest-ireighbor and next-nearest-neighbor 
interactions (J and J' respectively). The core physics 
can be generalized to the inclusion of further neighbor 
interactions. The bond ordered state is proposed to 
correspond to the spin gaped state which occurs when 
0.4 < J/ J' < 0.6 0, [13 • The columnar BD or- 
der renders the system to be equivalent to an effective 
coupled 2-leg spin ladder system (see Fig. 2). It is sup- 
posed that at moderate doping the BD and spin gap are 
retained P, which entails the confinement of spin-i 
excitations (i.e., pairing). The spin gap is the prominent 
energy scale which characterizes the BD. In the spirit of 
the RVB doped spin liquid ideaj^, pairing is inherited 
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FIG. 2: Columnar bond dimerization (BD); It is conjectured 
by Sachdev et al.^ 01 that the low energy physics is gov- 
erned by an effective Hamiltonian of a 2D t-J model where 
the Heisenberg interaction J is alternating in strength along 
one axis. 




from the spin liquid state of the undoped BD Heisenberg 
model. An analogy is drawn with the case of dimeriza- 
tion in one dimensional spin chains, where the spin gap 
can be as large as J (the Majumdar- Ghosh limit). A- 
priori, it might be that columnar BD can lead to a spin 
gap on the scale of weakly coupled 2-leg ladders ^ ^ . 

Since the undoped cuperates have a quasi long range 
antifcrromagnetic ground state 0, ^3 j the applicability 
of a BD model to the cuperates is not straight foreword. 
The mechanism implied by the phase diagram (Fig.l)j3| 
is the following. The microscopic degrees of freedom and 
interactions in the cuperate superconductors are capture 
by a one band generalized t-J model (e.g., t-t'-J-J'-J" ) 
plus coulomb interactions. The bare parameters values 
are such that groundstate of the undoped system is the 
Neel state as described by Chakravarty et al. 0. With 
the introduction of holes, the holes dynamics enhances 
the effective longer range antiferromagnetic interactions 
(i.e., imagine integrating out high energy hole hopping 
processes to arrive at an effective lower energy Hamilto- 
nian). Thus, at intermediate energy scales, it is conjec- 
tured that the effective spin interactions model parame- 
ters are in the fixed point regime of the bond dimerized 
Heisenberg model. The low energy holes dynamics is 
self consistently centered on the strong bonds (which in 
columnar BD form seem very much like 2-leg ladders) . 

The addition of coulomb interaction has the effect of 
changing the spacing of the hole rich columns compared 
with the model without such coulomb interaction 
Hence, a stripe like structure of alternating hole rich 
lines and hole poor domains can be formed (see insets in 
Fig.l). We comment that the BD as envisaged in [3, 01 is 
a property not only of the hole rich stripe but also of the 
hole poor environment between stripes. Thus, the spin 
gap in the hole rich stripe domain is tightly linked with 
properties of the hole poor domain between stripes. 



FIG. 3: Numerical simulation results for the spin gap of Ji — 
J2 Heisenberg models^^l- Size scaling of the energy gap to 
the first S = 1 spin excitation. Largest system is A'^ = 144. 
J2/J1 = 0.38 (triangles), 0.45 (squares), 0.50 (full circles). 
For comparison, the unfrustrated (J' = 0) Heisenberg model 
data is given (empty circles) . 



CRITICAL TEST OF THE PROSPECTS OF 
COLUMNAR DIMERIZATION 

Prior to discussing direct tests of the spontaneous BD 
mechanism, one may ask what are its most idealized 
prospects to begin with? Numerical simulations results 
of for the spin gap on a square Ji — J2 Heisenberg model 
are shown in Fig-3 10]. A spin gap is expected only in 
the range of parameters 0.4 < J2/J1 < 0.6. It appears 
that, unlike the one dimensional spin chain case where 
the spin gap can be as big as As ~ Ji (the Majumdar- 
Ghosh limit), the effective maximal dimerization in a 2D 
Ji — J2 Heisenberg model on a square lattice is never giv- 
ing rise to a gap value bigger than Ji/8. While it may be 
that by fine tuning some additional further interactions 
(>^37 Jii etc..) one might conjure higher gap values, that 
would not constitute a robust model. Since the square 
Ji — J2 Heisenberg model is the effective model system 
which is at the origin of the BD mechanism proposed 
by Sachdev, the results presented in Fig. 3 already put 
in serious doubt the basic conjecture which underlies the 
introduction of spontaneous BD as responsible for phe- 
nomena at energy scale A > J/6 in the cuperate high 
temperature superconductors. 

Even leg ladders are an ideal testing ground for t- 
J model physics. Numerical simulations of doped t-J 
model ladders 0, show groundstate pairing corre- 
lations in agreement with effective columnar BD (see 
Fig. 4). Therefore, ladders seem to be a favorable sys- 
tem for BD. It is worth noting that the appearance of 
columnar pairing correlations in ladders may also be a 
simple consequence of geometrical anisotropy and bound- 
ary conditions. Yet, such knowledge is not essential to 
our arguments in this paper. We argue that the key issue 
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FIG. 4: DMRG results for the real space singlet pairing expec- 
tation values on nearest-neighbor sites |l5ll . Pair correlations 
on doped 4-leg ladder seem to manifest columnar dimeriza- 
tion. 
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FIG. 5: DMRG results for the real space singlet pairing expec- 
tation values on nearest-neighbor sites flSt . Pair correlations 
on doped 4-leg ladder seem to manifest columnar dimeriza- 
tion. 



is the size of the spin gap. 

Numerical simulations by CORE method on doped 4- 
leg ladders already seem to indicate that the spin gap of 
the doped ladder does not get bigger than that of the 
undoped ladder (see Fig. 5). But there is an ambiguity 
between open and periodic boundary conditions. Appar- 
ently, the spin gap of 4-leg ladders (Aq « J/5) is to begin 
with too large to draw sharp conclusions from. On the 
other hand, as we elaborate below, similar calculations 
for 6-leg ladders are all that is needed to make the point. 

Undoped 6-leg t-J ladder models have a small spin gap 
Aq ~ J/20. Upon doping, the ground state still has a 
spin gap and enhanced pairing correlations. We propose 
to examine a 6-leg ladder model with the same bare t-J- 
J' model parameters as are assumed for the C-0 planes 
of the cuperate superconductors (e. g.. b are ring exchange 
interaction may need to be included The generation 

of effective interactions due to the holes dynamics is a 
local process which qualitatively and quantitatively we 
expect to be the same in wide ladders (e.g., 6-leg or 8- 
leg ladders) an in a 2D plane. Consequently, what ever 
is the potential for BD in the 2D C-0 plane, it should 
manifest itself with the same strength already in wide 
ladder models. 

The spin gap in undoped and doped 6-leg t — J ladders 



can be deduced from the energy gap to the first 5* = 1 
spin excitation. Examples of such a calculation were done 
for a 2-legand 4-leg ladders by both DMRG and CORE 
methods jl2tll^ . For spontaneous BD to be a mechanism 
relevant at energy scales of order A « 25QK, it is required 
that the size of the spin gap for the doped system will 
be on the order of Ag « 250K, i.e., 4 times bigger than 
that of the undoped 6-leg ladder!! 

In summary, bond dimerization (BD) is seen to occur, 
to some extent, in doped even-leg ladders. Thus, aspects 
of the resulting physics can be directly studied by nu- 
merical investigations of such ladders. Specifically, we 
have proposed testing whether the energy scales involved 
are sufficiently large to be relevant to a host of proper- 
ties of the cuperates (such as pairing and stripes). The 
experimental phenomenology is sufficiently complicated 
that detailed quantitative comparisons are certainly pre- 
mature, but we believe there is a significant theoretical 
distinction between the maximum spin-gap energy scales 
As w J/8 found in BD magnets 10], and the minimum 
necessary magnitudes A^ > J/6. More directly, if the 
BD mechanism fails to produce a high spin gap energy 
scale in doped 6-leg and 8-leg ladder systems (as is our 
prediction) then its prospects as the root cause of pairing, 
stripe formation, or normal state properties in the cu- 
perate superconductors are in serious doubt. Of course, 
we cannot rule out the possibility that additional well 
adjusted interactions, such as the electron-phonon inter- 
action, could increase the stability of the BD state and 
hence manage to overcome this objection. Yet, such ele- 
ments are not part of the framework currently presented 
byfliii. 

We thank Luca Capriotti, Doug Scallapino, and Steve 
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